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Abstract. - We propose a ballistic coalescence model (punctuated-Hamiltonian approach) 
mimicking the fusion of vortices in freely decaying two-dimensional turbulence. A temporal 
scaling behaviour is reached where the vortex density evolves like . A mean-field analytical 
argument yielding the approximation £ = 4/5 is shown to slightly overestimate the decay 
exponent £ whereas Molecular Dynamics simulations give £ = 0.71 ± 0.01, in agreement with 
recent laboratory experiments and simulations of Navier-Stokes equation. 



Two-dimensional turbulence has the fascinating property of organizing into coherent struc- 
tures from a disordered background. This feature has been observed in laboratory experiments 
[[[], |2| and in numerical simulations, as first emphasized and investigated by Mc Williams ||. 
It can account for the robustness of Jupiter's Great Red Spot, a large scale vortex which has 
been persisting in a turbulent shear for more than three centuries Q . 

We shall address the issue of freely decaying two-dimensional turbulence, for which the 
decay from random initial conditions can be divided into three stages. In the first stage, 
the system self-organizes into a set of coherent vortices containing most of the flow vorticity. 
The evolution is then dominated by the mutual advection of these structures, punctuated by 
dissipative events: whenever two like-sign vortices come closer than a critical distance, they 
merge (coalesce) to form a bigger vortex, while dipoles form a very stable state. Finally, the 
third stage starts when there are very few dipoles left, which decay diffusively. 

In this Letter, we shall concentrate on the second stage of the evolution, during which 
the average extension of the vortices and their relative distance grow whereas their number 
density n decreases in time. Numerical simulations of the Navier-Stokes equation || |(| have 
shown that this decay is algebraic: n(t) oc with £ ~ 0.72 ± 0.03. On the experimental 
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side, recent investigations in thin stratified layers of electrolyte emphasize the importance of 
coherent vortex dynamics for the decay of turbulence in two dimensions [0, || . 

Assuming a self-similar evolution of the vortex system, one can infer from dimensional 
grounds that n(t) oc t~ 2 , provided the energy is the only conserved quantity |^|, §]. However, 
the emergence of coherent structures significantly slows the density decay (£ < 2). In addition, 
numerical solutions J|, ^] indicate the appearance of a second conserved quantity: the vorticity 
amplitude lo inside the vortex cores. On the assumption that all the vorticity is concentrated 
in the vortices and given the two above mentioned invariants, the authors of J|, ^| derived 
a non-conservative scaling theory expressing all statistical properties in terms of the scaling 
exponent £, and constructed the simplest model capturing the essential features of the self- 
similar stage of freely decaying two-dimensional turbulence. As no inviscid invariants other 
than the energy and the typical vorticity are preserved, their theory is non-conservative. The 
vortices are modelized by discs of radii <7j, uniform vorticity u> and thus circulation I\ = it a 2 u. 
The equation of motion of vortex i, having position r,(t) = {xi(t), yi(t)}, is assumed to be 
governed by the Hamiltonian dynamics of point- vortices (the so-called Kirchoff 's laws |flO|): 



where the Hamiltonian Ti. takes the form 

H = \j2T i T j G{T i {t),T J (t)) (2) 

and G is the Green's function inverting the laplacian V 2 with the required boundary conditions 
{e.g. Dirichlet-like in a closed domain). The dynamics (preserving all moments of the vorticity) 
is punctuated by inelastic merging events with a rule derived from the conservation laws: 
since the energy £ of a vortex is of order T 2 oc uj 2 er 4 , the local conservation of £ with a time 
independent vorticity implies that the vortex resulting from the coalescence of vortices i and 
j has a radius <Xy given by 

a t] = (*l + *^ 4 . (3) 

This punctuated-Hamiltonian approach has numerically given a decay exponent £ ~ 0.72±0.02 
(cf |^|, p|), in agreement with the simulations of the full Navier-Stoke equation, indicating that 
the restriction of vortex interactions to those that have the most dramatic impact on the flow 
(i.e. those in which the number of vortices changes), seems sufficient to account for the decay 
exponent £. 

Several theoretical attempts have been made to determine £. Motivated by the scaling 
laws of H, Pomeau proposed an argument yielding £ = 1, but argued for lowering corrections 
fr| . In the possibly related context of Ginzburg-Landau vortex turbulence, the value £ = 3/4 
has been put forward On the other hand, mapping the original problem on a solvable 

model of charges on a lattice, Sire obtained n(t) oc ln(t)/i, where the logarithmic correction 
could explain why an effective exponent £ lower than 1 is measured [ p^[ . So far, the accuracy 
associated with system sizes available for the simulations or the experiments does not allow 
to discriminate such a law against a truly algebraic behaviour. 

In this Letter, we shall further simplify the original punctuated-Hamiltonian model of 
Carncvalc et ai, then obtain a simple analytical estimation for £ and test its validity by 
extensive Molecular Dynamics simulations. We modelize the motion of vortices between 
merging events by free flights, and concentrate on the asymptotic scaling regime followed by 
such a ballistic coalescence model, which will be shown to capture the essential features of the 
decay stage. A system initially made up of Nq identical discs of radius erg i n a two-dimensional 
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domain of surface S, with velocities Vj sampled from a distribution ipo(v), undergoes free 
particle motion, interrupted by a succession of instantaneous, completely inelastic collisions 
whenever a pair of discs comes into contact. As a result of a binary collision, the two initial 
discs merge into a single disc, which radius is given by eq (^). The corresponding merging rule 
for the velocity is needed to completely specify the model, but we shall indicate that its precise 
form does not affect the asymptotic properties of the scaling regime. A particular model of 
ballistic coalescence with conservation of mass and momentum has been introduced in juj . We 
first present a heuristic argument justifying the occurrence of an algebraic asymptotic regime 
on dimensional grounds. From the density n(t) = N(t)/S, the instantaneous average radius 
(a(t)}, the instantaneous average modulus of the velocity in the center-of-mass frame (v(t)}, 
and time t, we can form two dimensionlcss quantities and write 



T, 



f=HTJ w 

where / is an unknown scaling function and the Boltzmann mean collision time t b (t) is the 
individual mean time between two collisions suffered by any one particle: t b is the ratio of the 
time-dependent mean-free-path to (v) 

T ° = • (5) 

(a) I (v) is the time needed for a typical vortex to fly on the distance of a typical diameter, 
during which the system suffers no collision. As a consequence of eq. (||), the surface coverage 
(or packing fraction) is indeed a decreasing function of time, so that (a) eventually becomes 
much smaller than the mean-free-path. Correspondingly, on a time scale (a) / (v) <C T B (t), 
the evolution is free which indicates that t b is the only relevant time scale for the dynamics. 
Dimensional analysis then has to be reconsidered, and implies that t b scales like the physical 
time itself: t b oc t. 

During a time interval r B , each vortex typically collides, leading to a decrease of density 
which is of the order of the density itself: 

dt T B W 

Inserting t b oc t into (^|) , we obtain an algebraic time dependence for the density: this "critical" 
behaviour can be attributed to the uniqueness of a relevant time scale in this problem. We 
introduce two exponents £ and 7 to describe the asymptotic regime: n(t) cx t~^, (v(t)) oc V , 
the conservation of a 4 during encounters implying n(cr 4 ) cx t° and then (a) oc i^/ 4 . From the 
definition (|5|) and the relation t b oc t, one obtains the corresponding scaling relation: 

^ = 1 + 7- (7) 

A more rigorous derivation, based on a scaling analysis of the BBGKY-likc hierarchy governing 
the evolution of the distribution functions can be found in |l5). We emphasize that eq. (Q) 
is exact in the scaling regime and independent of the coalescence rule defining the velocity of 
a vortex resulting from a colliding pair. An other but approximate relation between £ and 
7 can be obtained from an analogy with a previously studied model of ballistic coalescence 
@ @. On that purpose, we introduce for each vortex a fake mass oc af and define a 
fake energy by rriivf. From eq. ([|), the mass is conserved during the collisions, which imposes 
that the average mass scales like the inverse of the density: (to) oc ifi. The total fake energy 
density n(mv 2 ) consequently scales like t~ 21 . The second and approximate relation between 
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£ and 7 is obtained assuming that each collision implies two typical uncorrelated particles, so 
that the fake energy loss equals the typical (i.e. average) fake energy. The fake energy per 
particle is then conserved ("isothermal" approximation), and the fake energy density has the 
same temporal decay rate as the number density: £ = —2j. Upon substitution into eq. (^), 
one obtains 

c = I- («) 

The above estimation may be coined "mean-field" and amounts to assuming that the linear 
momenta mjVj of the colliding vortices are uncorrelated p6| . 

Our estimate £ . = 0.8 is in fair agreement with the experimental and Navier-Stokes 
results, but needs to be confronted to numerical simulations. On that purpose, we specify 
a transformation rule for the velocity v„ of a vortex v resulting form the merging of vortices 
1 and 2, by coming back to Kirchoff's law. Equation ([[]) can be recast in the form 

v i = E r j v ^ -> !) + r 2 v ( 2 -> !) ( 9 ) 

v 2 - ^ r^v(j^2) + r lV (i->2) (io) 

where v(j — > i) is the velocity field created at point i by a vortex j of unit circulation. Just 
before the fusion, vortices 1 and 2 are close, and vortices j ^ 1, 2 are on average far apart form 
the pair 1-2. We thus consider 

v(i -» 1) ~ v(j -» 2) ~ v(j -> v) for 1,2. (12) 

Taking into account the symmetry constraint v(j — > fc) = — v(fc — > j), eqs. (|9|)-([TT1) supple- 
mented with approximation ( |l2| ) yield 

(ri + r 2 ) v„ = Tivi + r 2 v 2 v w = -g-^ — 5 (o-i vi + ^ v 2) • (13) 

+ 02 " ' 

We have performed two-dimensional Molecular Dynamics (MD) simulations of the ballistic 
coalescence model defined by merging rules (|) and (|l|), with an initial number of vortices 
5 10 3 < A^o < 510 5 . Periodic boundary conditions with a square simulation cell have been 
implemented and the equilibrium state of elastic hard discs served as the initial condition. After 
an initial transient corresponding to t < T B (t = 0), an algebraic scaling regime is evidenced, 
with a decay exponent £ ~ 0.71±0.01 (see figure 1) independent of iVo and of the initial packing 
fraction. The discrepancy between the mean-field prediction (£ . = 0.8) and the simulation 
result can be attributed to correlations in velocities required for the occurrence of a collision. 
These correlations are neglected in mean- field, which overestimates exponent 7, and from (^), 
overestimates £. Our exponent is very close to the Navier-Stockes value £ ~ 0.72. Moreover, 
the mean-free-path is found to scale like i 3 ^/ 4 , with a growth exponent 0.53, in agreement 
with the experimental value 0.45 ± 0.1 reported in M. These results indicate that the essential 
features of the decay stage can survive to a significant simplification of the complex interactions 
between vortices. 

The robustness of £ with respect to the merging rule for velocity has been tested by MD 
simulations with conservation of er 4 and of cr 4 v in collision events. This system is the analog of 
the original ballistic coalescence model introduced in Ej] , in the sense that the conserved linear 
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momentum <r 4 v involves a mass cr 4 which is itself a conserved quantity (this is not the case 
with rules and ( |T3| ) for which no conserved quantity involving the velocity can be defined). 
The model again displays an asymptotic density decay with £ ~ 0.70 ± 0.02, irrespective of 
the initial conditions. For the latter collision laws, the theory of Smoluchowski rate equations 
[ fl7| yields a mean-field approximation £ = 4/5 [ jl6| , confirming our estimate for £ mf . 

In summary, we have derived a ballistic coalescence model mimicking the fusion of vortices in 
freely decaying two-dimensional turbulence. The decay exponent £ obtained is in agreement 
with recent laboratory experiments |/| ^| and close to its Navier-Stockes counterpart. The 
precise form of the velocity merging rule is found to have little influence on £ which depends 
more crucially on the merging law for the size (eq. (|3|)). Our simplification of the original 
punctuated-Hamiltonian approach proposed by Carnevale et al. || further allows to put 
forward the simple mean-field analytical estimation £ . = 0.8, and indicates that conservation 
rule (0) is the crucial ingredient governing the self-similar stage of the decay. However, the 
evolution does not conserve the surface coverage 7m((7 2 )/4, which decreases like t~^ 2 : the 
mean vortex size eventually becomes much smaller than the typical distance nT 1 !" 1 (itself 
vanishingly small with respect to the mean- free-path) , and a vortex can interact with many 
neighbours between two collisions in such a way that its motion is likely to be quite erratic 
in the asymptotic regime. The validity of ballistic motion is thus questionable for very low 
packing fractions, whereas it seems more reasonable for high surface coverages (i.e. in the 
early stages of the dynamics). The correspondence between the value of £ obtained here and 
those reported previously ^ ^, fj], || could then point out to the non asymptotic character of 
the latter. Finally, unlike the original turbulence solution ||, our model does not conserve 
the typical velocity (v) (7^0 and we measure 7 ~ — 0.47± 0.01): we argue that the relevant 
velocity for our model is the relative velocity of close vortices (i.e. the velocity at the scale 
of a typical diameter), whereas the typical velocity of vortices associated with the conserved 
energy refers to the larger scale of inter-particle distances |l§| ]. As we expect an energy 
spectrum E(k) cx k~ a with a > 1, the ratio of these two quantities scales like 



( fc = V( g )) K rf(l-a)/8 



v(k = ^/n) 
and is thus a decreasing function of time. 



t«l-a)/8 (14) 
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10g(t/T„) 

Fig. 1. - logio(no/n — 1) vs. logio (£/ T o)> where to = T B (t = 0) and no is the initial surface density 
corresponding to a packing fraction 0.24. The initial number of vortices is iVo = 10 5 . The dashed and 
dotted lines have respectively slopes 0.8 (mean-field prediction) and 0.71 
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